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Abstract. We describe our first attempt at modelling nucleosynthesis in mas- 
sive AGB stars which have undergone core carbon burning, the super- AGB stars. 
We fit a synthetic model to detailed stellar evolution models in the mass range 
9 < M/Mq < 11.5 (Z = 0.02), and extrapolate these fits to the end of the 
AGB. We determine the number of thermal pulses and AGB lifetime as a function 
of mass and mass-loss prescription. Our preliminary nucleosynthesis calculations 
show that, for a reasonable mass-loss rate, the effect of hot-bottom burning in 
super- AGB stars on the integrated yield of a stellar population is not large. There 
are many uncertainties, such as mass-loss and convective overshooting, which pre- 
vent accurate yield calculations. However, as potential progenitors of electron- 
capture supernovae, these stars may contribute 7% of non-type-la supernovae. 
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Introduction 

Stars are traditionally divided into those 
which explode and those which do not. 
Stars which explode as supernovae by a 
core-collapse mechanism are the massive 
stars, while their lower-mass cousins en- 
ter a thermally pulsing asymptotic giant 
branch (TPAGB) phase where rapid mass- 
loss ends their evolution. It is difficult to 
specify the exact mass boundary between 
these popu lations, with estima tes ranging 
from 7 M Q iGirardi et all2Q0Cl) l in to more 
than 11 M ? l|Ritossa et alJIlffflfl . depend- 
ing on metallicity and convective over- 
shooting. 



Send offprint 

R . G . IzzardOphy s .uu.nl 



requests 



to: 



The problem relates to the fate of the 
star after carbon ignition, which depends 
on the degeneracy of the core. In high-mass 
stars the non-degenerate core burns car- 
bon, then neon, oxygen and silicon. A core- 
collapse supernova soon follows, leaving a 
neutron star or black hole. A lower-mass, 
but massive enough to ignite carbon, (par- 
tially) degenerate co re burns ca rbon in a 
series of flashes (e.g. lSiesdl2006j) and then 
moves to a super thermally pulsing AGB 
(SAGB) stage, with double-shell burn- 
ing above a degenerate oxygen-neon core, 
where mass-loss terminates the evolution, 
leaving a white dwarf (Siess, this volume). 
In a small number of stars the oxygen- 
neon core may grow beyond 1.368 Mq dur- 
ing the pulsing phase, which leads to a col- 
lapse of the core due to electron capture on 
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24 Mg (Poelarends, this volume). Detailed 
models of these stars have been con- 



structed by Garcia-Berro k Ibenl ll 19941) 



iRitossa et al.l lll996h: l(^arcia^Berro et al.l 
lib en et all (^9J|) and more re- 



cently bv lEldridee k Tout! l)2004|) who con- 
sider SAGB stars as observable supernova 
progenitors (see also Poelarends, this vol- 
ume). 

Once carbon ignition and second 
dredge-up have finished, SAGB stars pulse 
in much the same way as normal TPAGB 
stars, albeit with a shorter interpulsc 
period. By virtue of their high mass, 
one expects high temperatures at the 
base of the convective envelope and as- 
sociated hot-bot t om b urning (HBB; e.g. 
iBoothrovd et alJ ^995). Their envelopes 
should be processed by the hydrogen- 
burning CNO, NeNa and MgAl cycles. 
For studies of Galactic chemical evolution, 
these stars may be an important source of 
nitrogen, sodium and aluminium and they 
may play a part in the globular cluster 
Na — O anti correlation mystery (D'Anto na, 
this volume: fVentura fc D'Antonall2005|) . It 
is unclear whether third dredge-up occurs 
in SAGB stars - current models suggest ei- 
ther it does (Doherty, this volume), does 
not (our mode ls and Siess, this v olume), 
or is inefficient IRitossa et al .l l1 996ft . There 
are no detailed studies of chemical yields 
from SAGB stars, probably because de- 
tailed stellar models take a long time and 
are difficult to construct, and suffer from 
the usual uncertainty due to mass-loss and 
convective overshooting. A synthetic mod- 
elling technique speeds up modelling and 
enables us to explore the uncertain param- 
eter space. 

In this paper we calculate the chem- 
ical yields of SAGB stars using a syn- 
thetic model based on the AGB model of 
llzzard et all ((20041 104). We approximate 
stellar structural variables with formulae 
and interpolation tables, and use a simple 
model for HBB to follow the CNO, NeNa 
and MgAl cycles and surface abundances 
of C, N, O, Ne, Na, Mg, Al and Si. The 
synthetic model is then used to extrapolate 



evolution beyond the detailed models to the 
end of the SAGB phase. This is possible 
because the structure of AGB stars is such 
that after a number of pulses, the evolutio n 
reaches a limit cycle l)Ritossa et alJ ll996). 
We calculate supernova rates and chemical 
yields, and also change the input physics 
(especially the mass-loss rate) to determine 
the effect of uncertainties. Such a parame- 
ter space exploration is currently impossi- 
ble with a normal stellar evolution code, 
because the CPU time required is simply 
too large. 

Models, full and synthetic 

Our full evolution models were con- 
structed with the STERN code 
(|Heger. Langer. k WoosTevI l2000|) . We 
constructed models of mass 8.5, 9.0, 10.0 
and 11.5 M Q which undergo 12,26,10,16 
pulses respectively, with metallicity 
Z = 0.02, no convective overshooting, 
no mass loss and no rotation. Further 
description of these models can be found 
in Poelarends et al. (this volume). 

Our synthetic models are based on 
those of 104 with some updates and changes 
for the SAGB phase (a detailed descrip- 
tion will be found in a later paper; Izzard 
and Poelarends, in preparation). The lumi- 
nosity formula was altered to fit our de- 
tailed models, and the radius follows from 
logi? ~ logL. The initial and post-second- 
dredge-up abundances, core mass at the 
end of core helium burning and core mass 
at the first thermal pulse, are interpolated 
from tables based on the detailed models. 
Stars with a helium core mass above 1.6 M Q 
during the early AGB ignite carbon, while 
stars with a degenerate oxygen-neon core 
with mass greater than 1.38 M© collapse to 
neutron stars. We assume there is no third 
dredge up. 

The synthetic HBB model was de- 
scribed in 104 but we use the latest version 
which includes the NeNa and MgAl cycles 
as well as CNO. It approximates the burn- 
mix-burn-mix. . . cycle in a real convective 
envelope (which has a thin HBB shell at 
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the base) as a single burn-mix event during 
each interpulse period, with a large fraction 1 
of the envelope burned for a given time. 
The temperature and density at the base 
of the envelope are fitted to simple formu- 
lae of the form 1 — exp (—Ntp) which ap- 
proaches a constant as Ntp , the number of 
thermal pulses, increases. The fraction and 1 
the burn time are calibrated as a function 
of stellar mass, as in 104. 

We apply one of the following mass- 
loss prescriptio ns during the SAGB phase : 
none, original IVassiliadis fc Woodl (1993, 
VW93), VW93 but lKa,ra,ka,s et all (|2002i 
K02) variant, R ei mers with 77 = 1 or 77 = 5 
llReimerslll975|) or iBlocker fc Schonbernerl 




1991) with 77 = 0.1. Prior to the SAGB, 
we either ap ply no mass loss, o r the 
compilation o flHurlev. Tout, fc Polsl <|2002t 
H02). In the cases where mass loss does 
not expose the core before it grows to 
1.368 Mg, the core is quietly converted 
into a neutron star while the envelope 
is ejected to space (this is the sub- 
ject of some debate e.g. iNomotol 119871 
(Gutierrez. Canal, fc Garcia-Berroll2005j) . 

Results 

We are confident that our extrapolation 
of the stellar structure (luminosity, radius, 
mass and core mass evolution) is reason- 
able, within the uncertainty that is mass 
loss, because a very similar model works 
well for lower-mass AGB stars (104). The 
mass-loss rate for SAGB stars is unknown, 
so we consider all the possibilities. Figure 
^ shows the number of thermal pulses dur- 
ing the SAGB phase which ranges from 30 
(Blocker, with SAGB lifetimes of 4, 000 to 
5 x 10 4 years) to 4,000 (Reimers 77 = 1, 
with lifetimes of 10 5 to 10 6 years), or 8, 000 
with no mass loss. Mass loss during core 
helium burning, prior to the SAGB, is in- 
cluded in the H02 prescription and is not 
negligible for stars above 8 M - they lose 
around 0.5 M Q during this phase, which af- 
fects the subsequent evolution and HBB. 

One way of constraining the value of 
mass- loss rate is through supernova counts. 



Fig. 1. Number of thermal pulses as a 
function of initial stellar mass for differ- 
ent mass-loss prescrip tions: no wind at al l 
or pre-SAGB wind of iHurlev et"afl l|2002j) 
with Reimers (77 = 1 or 5), VW93 or 
Blocker mass-loss during the SAGB. 



If there is little mass loss during the SAGB, 
many stars' cores should reach the elec- 
tron capture limit of 1.368 M Q before their 
envelope is lost. In the limit of no mass 
loss during the SAGB, we find the ratio of 
the electron-capture to type-II supernova 
rates is about one 1 . With the VW93 (K02 
variant) mass-loss the ratio drops to 7%. 
If electron capture supernovae are distin- 
guishable from normal core-collapse types 
(and accretion induced collap ses in bina- 
ries ;|Podsiadjomk^^3in3) then we can 
constrain the mass- loss rate 2 . 

In figure El we show the results of our 
HBB calibration and extrapolation for a 
10 Mq star. The detailed models truncate 
their abundance output to the nearest 10~ 3 
in the log, which causes the stepping be- 
haviour, while the stepping in the synthetic 
model is because HBB is done at the end 
of each pulse. Ten pulses is just enough to 
calibrate the isotopes which change rapidly, 
13 C, 17 and 21 Ne are also useful in this re- 



1 We do not include an upper mass limit for 
SNe, so ratios quoted are lower limits. 

2 We only model single stars here, binary in- 
teractions will alter the result. 
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Fig. 2. Log surface abundance by mass fraction vs time for 12 C and 14 N for our 10 M Q 
detailed (black lines) and synthetic (grey lines) models, both with no mass loss. The 
top panels show the result of our HBB calibration, the bottom panels the same models 
extrapolated to the end of the SAGB. 



gard, while other isotopes such as 16 O are 
not. We have to assume the other species 
( 20 > 22 Ne, 23 Na, Mg and Al) follow from this 
calibration. The situation is little better in 
the 9 Mq star: it has 26 pulses, but its lower 
temperature slows the burning. 

Finally, we consider the magnesium iso- 
topes. Figure OH shows the surface abun- 
dance of 24 Mg and 25 Mg as a function of 
time for 9 < M/M Q < 12 with H02 mass 
loss prior to the SAGB, and VW93 (K02 
variant) mass loss during the SAGB. Only 
the 11 and 12 M Q models show significant 
burning of 24 Mg to 25 Mg, with up to a fac- 
tor of two increase. For M < 10 M Q the 



HBB is simply not hot enough to enable 
the MgAl cycle. At all masses, rapid mass 
loss turns off the MgAl cycle after about 
2 x 10 4 years 3 . 

Discussion 

This was our first attempt to model these 
stars synthetically and, at least for the 
structure variables such as luminosity, ra- 
dius, core mass (and growth), we have con- 
fidence in our model. However our HBB 

3 As in 104 we modulate the temperature 
with a factor (JV/env/iWenv,iTp) ' 02 to turn off 
HBB as mass is lost. 
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Fig. 3. The surface abundance of 24 Mg and 25 Mg during the super- AGB phase for the 
initial mass ranee 9 < M/Mq < 12. Hot-bottom burning converts 24 Mg into 25 Mg, while 
mass-loss stops the burning and terminates the evolution. 



model is not as solid, given that we have 
extrapolated forward by a factor of ten 
or more in time compared to our de- 
tailed models. It is clear that we must ex- 
tend our detailed models to at least cover, 
say, half the evolution (which may still be 
many thousand pulses), and see if the syn- 
thetic model predictions match the detailed 
model 4 . 

Even when we have extended the de- 
tailed models, we will still suffer from the 
mass-loss uncertainty. It seems, from fig- 
ure m we can get as many pulses as we 
like. New observations are providing in- 
sight into the mass-loss rates in oxygen-rich 
AGB stars and we integrate these new rates 
into o ur next calculations ijvan Loon et al.l 
2005). We have also neglected the problem 
of convective overshooting, which reduces 
the mass for formation of SAGB stars by 
about 2M© (Siess, this volume). The ini- 
tial mass function is a steep power law in 
mass, so a lower mass limit means more 
SAGB stars, and their lifetimes will be 
longer (lower mass means smaller luminos- 
ity, radius and M). 

Our preliminary results suggest that 
SAGB stars are not very important ei- 

4 The numerical problems seen in figure [5] 
will become irrelevant in that case. 



ther for Galactic or globular cluster chem- 
ical evolution. The bulk of element pro- 
duction comes from either lower-mass AGB 
stars or higher-mass stars and type II su- 
pernovae. This conclusion may change if 
SAGB stars suffer either third dredge-up, 
or dredge-out, a phenomenon where the 
carbon flash causes helium-b urned material 
to be mixed to th e surface l)Ritossa et al.l 
119991: ISiessi H)06). Also, we only con- 
sider envelope ejection for a collapsing 
oxygen-neon core. In reality, some of the 
core may be ejected too, particularly if 
some carbon remains after the core flash es 
l|Gutierrez. Canal, fc Garria-BerrcfeOQSD. 

Conclusions 

In the context of galactic chemical evo- 
lution, chemical yields from hot-bottom 
burning SAGB stars are not important, at 
least for most isotopes. If mass-loss rates 
are much lower than those expected from 
extrapolation of normal AGB rates, if there 
is dredge-up or dredge-out, or if our simple 
extrapolations fail, this conclusion may be 
premature. We are working to extend our 
detailed model set to remove the extrapo- 
lation problem. 

The ratio of electron-capture to type- 
II supernovae in single stars is about 7% 
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if we assume a IVassiliadis fc Woodl l|l993l) 
wind for the SAGB phase - this should be 
detectable, if it is possible to distinguish 
electron-capture from core collapse super- 
novae. 
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